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5 PRIORITY INFORMATION 

This application claims priority from provisional application Ser. No. 60/413,799 
filed September 26, 2002, which is incorporated herein by reference in its entirety. 

BACKGROUND OF THE INVENTION 

10 The invention relates to the field of Cherenkov radiation, and in particular to 

using photonic crystals as a medium to exhibit anomalous Cherenkov radiation. 

Cherenkov radiation (CR) is the coherent electromagnetic response of a medium 
driven by the swift passage of a charged particle. It is thus an effect strongly dependent 
on the medium dispersion. In a uniform, isotropic medium with frequency-independent 

15 permittivity e and permeability |i, the condition for CR is well-known, where the velocity 

of the particle v must exceed the phase velocity of the medium u ph = cj yfs/u . For a 
dispersive medium, such as a nonmagnetic material with a Lorentz-form dielectric 
response s(co), u P h(co)= cj yje{co) , is a function of frequency co. Because s(co) can 

reach arbitrarily high values near a resonance, it was long recognized that CR in a 1 
20 dispersive medium can happen for small charge velocities, e.g., o<u ph (0). The sub- 
threshold CR in a material near its phonon-polariton resonance was disclosed in 
Cherenkov Radiation at Speeds Below the Light Threshold: Phonon-Assisted Phase 
Matching, by T.E. Stevens, J.K. Wahlstrand, J. Kuhl, and R. Merlin, SCIENCE, vol. 
291, No. 5504 (26 Jan. 2001). 
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On the other hand, it was conjectured that there could exist another class of 
materials which have both 8 and jj. being negative, henceforth referred to as "negative 
index materials." The properties of such materials were disclosed in The 
Electrodynamics of Substances with Simultaneously Negative Values of s and ju , by 
V.G. Veselago, SOVIET PHYSICS USPEKHI, vol. 10, No. 4 (Jan. - Feb. 1968). It was 
suggested that a negative index material would reverse many of the well-known laws of 
optics. In particular, CR effect is predicted to be reversed, i.e., a fast-moving charge in a 
negative index medium should radiate in the direction opposite to that of its velocity. 

A further possibility exists when the charged particle travels near a periodic 
structure, where simple Bragg scattering can give rise to radiation without any velocity 
threshold. This phenomenon (the Smith-Purcell effect) was disclosed in Visible Light 
from Localized Surface Charges Moving across a Grating, by SJ. Smith and P.M. 
Purcell, PHYSICAL REVIEW, vol. 92, No. 4 (15 Nov. 1953). The radiation due to 
traveling charged particles has since been studied in one-dimensionally periodic 
multilayer stacks in Cerenkov Radiation in Inhomogeneous Periodic Media, by K.F. 
Casey, C. Yeh, and Z.A. Kaprielian, PHYSICAL REVIEW, vol. 140, No. 3B (8 Nov. 
1965), and near the surface of dielectric structures in Interactions of Radiation and Fast 
Electrons with Clusters of Dielectrics: A Multiple Scattering Approach, by FJ. Garcia 
de Abajo, PHYSICAL REVIEW LETTERS, vol. 82, No. 13 (29 March 1999). 



SUMMARY OF THE INVENTION 

According to one aspect of the invention, there is provided a system for exhibiting 
Cherenkov radiation. The system includes a beam of traveling charged particles. A 



photonic crystal structure receives the beam of charged particles. The charged particles 
move in the photonic crystal structure so that CR is produced at all velocities without 
requiring resonances in the effective material constants of said photonic crystal structure. 
According to another aspect of the invention, there is provided a method of 
5 exhibiting Cherenkov radiation. The method includes providing a beam of traveling 
charged particles. Also, the method includes providing a photonic crystal structure that 
receives the beam of charged particles. The charged particles move in the photonic 
crystal structure so that Cherenkov radiation is produced at all velocities without 
requiring resonances in tKe effective material constants of said photonic crystal structure. 

10 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic diagram of a radiation mode in a photonic crystal; 
FIG. 2 is a schematic diagram of a radiation cone; 

FIG. 3 A is a TE bandstructure of a 2D square lattice with air columns; FIG. 3B is 
15 schematic diagram illustrating a moving charge in the photonic crystal; 

FIGs. 4A-4C are graphs illustrating the relationship between the frequency, the 
wavevector k, and the corresponding group velocity u of the radiation modes in the first 
band of the photonic crystal of FIG 3B; 

FIG. 5A-5D are schematic diagrams illustrating the results of the radiation of a 
20 charge moving in the photonic crystal of FIG. 3B with v= 0.1c, 0.15c, 0.30c, and 0.6c; 

FIGs. 6A-6D are graphs showing the distribution of the radiation magnetic field 
perpendicular to the 2D plane with v= 0.1c, 0.15c, 0.30c, and 0.6c; 
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FIGs. 7A-7D are graphs showing the frequency spectrums of the radiation flux 
along the z-axis with v= 0.1c, 0.15c, 0.30c, and 0.6c; 

FIG. 8 is a schematic diagram illustrating a particle detector; and 
FIG. 9 is a schematic diagram illustrating a radiation source. 

5 

DETAILED DESCRIPTION OF THE INVENTION 

The invention describes using photonic crystals, i.e. periodic lattices for 
electromagnetic waves, as a medium which exhibits Cherenkov radiation (CR) and 
proposes new devices that make use of this effect. When a photonic crystal is coherently 

10 excited with a beam of traveling charged particles, the CR is produced. The radiation has 
its origin in both the transition radiation, which occurs when the charge experiences an 
inhomogeneous dielectric environment provided by the crystal, and the conventional CR 
in a uniform material, in which coherence is preserved throughout the medium. Unlike 
the Smith-Purcell effect, in which light is generated near a periodic grating and then 

15 propagates down through a uniform medium, the CR in this invention is generated and 
propagates through the same crystal in the form of Bloch waves. 

Due to the very complex Bragg scattering effect and the rich photon dispersion 
relations in a photonic crystal, the properties of Bloch waves can be very different, from 
waves in a uniform medium, leading to a variety of unusual phenomena. In particular, a 

20 charge moving in a photonic crystal radiates at all velocities without requiring the usual 
Cherenkov threshold condition or resonances in the effective material constants. 
Moreover, the invention predicts new CR wavefront patterns that are impossible to 
achieve within either a uniform medium or the Smith-Purcell effect. Furthermore, the 
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invention demonstrates situations in which CR propagates backward, a behavior 
reminiscent of that predicted in negative-index materials, without requiring the material 
constants to become negative. The invention provides systematic methods for analyzing 
these new phenomena, and confirms predictions by direct numerical simulations. 
5 By discussing the general condition for Cherenkov radiation in photonic crystals, 

it is important to consider a particle of charge q moving at a constant velocity v on a z 
axis inside a photonic crystal. An analytical expression for the fields generated by such a 
moving charge can be derived using the standard normal-mode expansion in Fourier 
space as: 

4^((k + g>v)gY«M-,) - i(k+g) .v, E l 

k«g ^(^-((k + g^v) J 

Here, the wavevector k is summed over all points in the first Brillouin zone, the band 
index n runs over all bands, and the reciprocal-lattice vector g is summed over all the 
reciprocal space. = £ G e knG e' (K+G) "' is the eigenmode of the photonic crystal at Bloch 
wavevector k and band index n with corresponding Fourier component e teG and eigen- 

15 frequency co^ (the photonic band structure), and is normalized according to 
l/Qjdr si^E^E'^ = S(k-l)S nm with Q being the volume of a spatial unit cell, of the 
photonic crystal. N is the number of the unit cells of the crystal, and NQ is the total 
volume under consideration. 

In the well-studied case of a uniform medium in which simple analytical 

20 expressions are available for £ ta and co^ , Eq. 1 can be shown to be equivalent to the 
familiar Fermi results. In the case of a general photonic crystal, the far-field radiation 
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modes can be deduced from the poles in Eq. 1. If restricted to positive frequencies, one 
can find the condition for radiation to be co^ = (k + g). v . Each solution set of k, n and g 
then corresponds to a CR mode emitted by the photonic crystal. The relative excitation 
strength of each mode is proportional to the magnitude of e te( } multiplied by slow 

5 functions of k and g on the numerator of Eq. 1. Since co^ = <0 (k+g)/J , the radiation 
condition may also be written as 

co^ = k»v Eq. 2 

with k now in an arbitrary "Briflouin zone, in- this way, g is just the reciprocal-lattice 
vector required to reduce k to the first Brillouin zone. 

10 Eq. 2 is the general condition for CR in a photonic crystal, and can be easily 

reduced to the usual threshold condition mentioned earlier in the case of a uniform 
dielectric medium. In a general photonic crystal, it may be solved by intersecting the 
plane k • v = co in k-space with the photonic-crystal dispersion surface co^ = co in the 
periodic zone scheme, as shown in FIG. 1. The dispersion surface is similar to a normal 

15 surface in crystal optics, but is the surface formed by k instead of the phase velocity. For 
each fixed n, solutions to Eq. 2 exist for arbitrary v. This is because, as co goes from the 
minimum to the maximum frequency of band n, the dispersion surface traverses through 
all possible k-points, and the plane k z =co/u only sweeps through a region with finite 
width in k z in k-space. Thus, no threshold exists for CR in photonic crystals. In 

20 particular, a moving charge can radiate even for v -> 0 . Physically, this effect is similar 
to the Umclapp process in solid-state physics. As u-»0 the solution for smallest 
frequency exists at co « g 0 *v , g 0 being a primitive reciprocal-lattice vector with smallest 
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positive g 0 «v . Note that since the radiation frequency goes to zero as u 0 , the present 
effect is fundamentally different from a similar result in dispersive medium, which 
requires e -» oo and radiation can only be near a fixed resonance frequency. In principle, 
this applies to any g, but the coupling strength, which is proportional to |e M _ g) |, must go 

5 to zero for large g. 

Positive frequency, as described in Eq. 2, requires that k^v > 0. In a uniform 
medium with positive index, energy flows along k, and thus a moving charge emit energy 
forward, i.e., the direction of radiation makes an acute angle with v. In a negative index 
material, the situation is reversed. The angle between energy flow and the velocity is 
10 now obstuse, and a moving charge must radiate backwards. 

In the case of a general photonic crystal, the direction of energy flow should be 
determined from the direction of the group velocity, i.e. the Poyning vector associated 
with the solution to Eq. 2. In the periodic zone scheme, there are many bands which have 
group velocities opposite to that of k. If these bands are coupled to an appropriate v in 
15 Eq. 2, then there is a CR in the direction opposite to that of v. Therefore, a negative 
group-velocity is a sufficient condition for reversed the CR effect. Note that this 
condition can be applied in the long-wavelength limit in the second Brillouin zone, where 
the photonic crystal exhibits positive effective index. 

Another characteristic of CR is the shape of the radiation wavefront. In a uniform 
20 nondispersive medium, it is well-known that a forward-pointing shock-front is produced 
on a cone behind the charged particle (the presence of dispersion removes the 
singularity). The difference between group velocities and phase velocities also alters the 
shape of the radiation cone. Here the invention presents a graphical method for 



determining the shape of the cone, as shown in FIG. 2, which can be easily extended to 
the case of photonic crystals. FIG. 2 illustrates a plot having both v and the group 
velocities u of all modes in Eq. 2 in a velocity-space. The trace of u is analogous to the 
ray surface in crystal optics, though in the present case it is formed by the group 

5 velocities of different frequencies. The magnitude of u is proportional to the distance 
traveled by the wavefront of the corresponding mode, and the magnitude of v is 
proportional to distance traveled by the moving charge. Thus, for each mode, the 
wavefront of radiation still lies on a cone whose angle is determined by the angle 6 
between v - u and v. The^maximum of all 0 f s determine the angle of the overall cone 

10 that encompasses all radiation. Propagating Bloch modes can only exist on the rear side 
of this overall cone, whereas the fields are evanescent on the forward side, and across the 
cone the radiation field amplitude experiences a drop. 

From these considerations it is understood that the group velocities of the modes 
determined by Eq. 2 are a key to understanding CR in photonic crystals. They can be 

15 calculated either from the gradient vectors to the dispersion surface or from the 
Heynman-Feynman Theorem specialized to photonic crystals. The detailed behavior of 
group velocities in a specific photonic crystal is studied further hereinafter. 

For concreteness, a specific problem in 2D is addressed to illustrate the- above 
analysis. Suppose a charge moves inside an infinite 2D square lattice of air columns in 

20 dielectric s = 12 , as shown in FIG. 3 A. The columns have a radius of 0.4a, where a is 
the period of the lattice, and the charge is moving in the 2D plane in the (10) lattice 
direction. The photonic-crystal TE (in-plane electric field E, appropriate for CR) 
bandstructure and the geometry of this problem are shown in FIG. 3B. Eq. 2 is solved to 
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determine the radiation modes of this photonic crystal and their group velocities, using 
software that calculates the photonic bands by preconditioned conjugate-gradient 
minimization of the block Rayleigh quotient in a planewave basis. For purposes of 
discussion, only the frequencies in the first band are analyzed and those in higher bands 
5 can be discussed similarly. The results are plotted in FIGs. 4A-4C for several different 
Vs. 

It is observed that the phase velocity of this photonic crystal in the long- 
wavelength limit (&>->0) is u c = 0.44c. For v«vc, the radiation coalesces into 

resonances around a> « G • v . For larger v, the resonances merge together to form 
10 emission bands outside which CR is inhibited (FIG 4A). As v increases, k and u within 
each emission band are strongly influenced by the photonic band structure (FIGs. 4B- 
4C). In particular, for v < v c most solutions exhibit negative z-components of group 

velocities, indicating backward-propagating Cherenkov effect for v < u c in general. 

Since radiation is allowed even for v «: o c , some modes can have group velocities that 

15 exceed v . According to the analysis of FIG. 2, modes traveling faster than the moving 
charge produce either nearly isotropic wavefront in all directions ( v < 0.1c), or reversed 
radiation cones that point backward (v = 0.15c). While the Smith-Purcell radiation 
always has a similar isotropic wavefront, it cannot produce a reversed radiation cone. 
The reversed cone effect is also forbidden by causality in any uniform passive medium, 
20 and is therefore a feature unique in photonic crystals. 

To confirm these predictions finite-difference time-domain simulations are 
. performed of the CR in photonic crystals, using a thick boundary perfectly matched layer 
(PML) which consists of 10 periods of identical photonic crystal buried inside. This 
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boundary condition can effectively absorb radiation modes in photonic crystals away 
from the band edge. The moving charge is implemented as a dipole of constant 
amplitude which points toward the direction of motion z and whose position depends on 
time. To verify negative CR, the fields are recorded on all points of an observation line 
5 perpendicular to z and the frequency components of the z-flux are calculated from those 
of the observed field values. 

FIGs. 5A-5D show the calculated group velocities and expected radiation cone 
shapes. The simulated radiation field pattern and the frequency spectrum of the z-flux 
are shown in FIGs. 6A-6D and 7A-7D, respectively. FIGs. 6A-6C demonstrate the 
10 radiation below u c , and FIGs. 7A-7C exhibit the negative z-flux frequency regions which 

are in agreement from what can be obtained from Eq. 2. Moreover, FIG. 6A shows the 
near-isotropic radiation at v = 0.1c , and FIG. 6B confirms the reversed CR cone. From 
FIG 6, the angular distribution of radiation is also readily available. When v=0.1c or 
0.1 5c, the radiation is distributed over a wide range of emission angles without producing 
15 a cone of intensity maxima. For v^O.Sc or 0.6c, however, the CR becomes collimated, 
and a definite emission angle in both the forward and the backward direction for most of 
the radiation energy can be observed. Note in FIGs. 6D and 7D that the CR radiation at 
high velocities (v>u c ) in the first band are qualitatively similar to that in a uniform 

medium with very small dispersion (i.e. forward-pointing radiation cone and positive z- 
20 flux). 

An important difference between radiation patterns for v <o c and v > u c is that 

the DC components (a> = 0 ) extends beyond the radiation cone in FIGs. 6A-6C, while in 
FIG. 6D the field outside the radiation cone is strictly zero. These DC components are 
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also reflected in the strong peaks around co = 0 shown in FIGs. 7A-7C. The finite flux 
values at co * 0 before the onset of negative z-flux is an artifact of finite computational 
size and time, and the DC components should not influence regions of co * 0 in ideally 
finite systems. Finally, it can be noted that there are strong, high-frequency "tails" of 
5 radiation behind the charge in all cases shown in FIGs. 6A-6D. These correspond to the 
radiation in the higher bands which have much smaller group velocity. Of course, here 
these tails can exhibit a backward radiation effect as well. 

Although a 2D "air-holes-in-dielectric" structure has been previously considered, 
the same physics obviously applies to 3D photonic crystals in general with little change. 

10 In particular, dielectric cylinders or spheres in air appear to be a good candidate for 
experimental studies of these effects. The charged particles can be chosen to be 
electrons, which with high velocities should be available from emission through a voltage 
difference of several tens or hundreds of kilo-volts. They can then be directed to travel in 
the all-air channels of these photonic crystals. Direct experimental verification of the 

15 anomalous CR effects in photonic crystals introduced in this invention should thus be 
very possible. 

Particles traveling at speeds below the phase-velocity threshold cannot be 
detected by conventional CR counters, and currently their observation relies on other 
devices, such as scintillation counters, proportional counters, or cloud chambers. These 
20 other devices, however, lack the unique advantages of strong velocity sensitivity and 
good radiation directionality as in conventional CR. With a photonic crystal, one should 
be able to achieve velocity selectivity and distinctive radiation patterns without any 
velocity threshold. Using the invention, one can form a particle detector 20, as shown in 
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FIG. 8, for counting of charged particles can now be made at arbitrary velocities using 
sensitive radiation detection devices, such as photomultipliers or vacuum photocell 
coupled to an amplifier. Charge particles 26 of all velocities are received in the all-air 
channel of a photonic crystal 28. The CR emitted by the charges is received by the 
5 photomultipliers 22 surrounding the photonic crystal 28, and the signals 30 can be 
amplified in an external amplifier. Note that waveguides may also be integrated into the 
photonic crystal to facilitate radiation collection. The sensitivity of this detector 20 only 
relies on the sensitivity of the photomultipliers 22 and does not require a velocity 
threshold. A number of devices measuring particle number, speed, charges or density 

10 flux can be built based on the detected radiation spectrum and wavefront pattern. Such 
devices are useful wherever particles of arbitrary energy are produced, such as in 
monitoring and controlling of nuclear reactors. 

Moreover, a radiation source design 36 can be formed using CR, as shown in FIG. 
9. An intense electron beam 44 is generated from high-voltage cathode 40 emission and 

15 pass through the all-air channel of a photonic crystal 42. It radiates and the radiation 46 
can be collected by waveguides 48 that are built inside the photonic crystal 42. Of 
course, the radiation can also travel as bulk Bloch waves as discussed above and be 
collected by external devices without using output waveguides. The frequency * of the 
radiation depends on its velocity and is tunable through the potential applied to the 

20 electron beam 44. Electromagnetic wave of arbitrary frequencies can now be generated 
using the effects as shown in FIG. 9. The range of the frequency of radiation 46 is now 
set by the photonic crystal 42 and the transparent spectral regime of the dielectric, thus 
selectively scalable beyond the optical wavelengths. Moreover, the frequency is tunable 
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at will by changing the velocity of the particles. Furthermore, integrated photonic crystal 
line-defects can be employed to concentrate and guide the radiation. A "dense" electron 
beam may also lead to coherence in the output radiation. This type of radiation sources 
allows very flexible designs, and will be especially attractive for frequencies that are 
otherwise difficult to access, e.g. in terahertz regime. 

All conventional applications of CR in the high- velocity regime, such as velocity 
measurement and selection, should be able to benefit from this invention. Since the 
particle can be made to travel in the air channel, the disturbance to the particle motion 
due to impurity scattering andrandom ionization, which introduce inherent losses in a 
conventional high-index material for CR detectors, is now completely absent. This could 
greatly enhance the performance of present CR detectors. 

Although the present invention has been shown and described with respect to 
several preferred embodiments thereof, various changes, omissions and additions to the 
form and detail thereof, may be made therein, without departing from the spirit and scope 
of this invention. 

What is claimed is: 
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